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^ We investigate nonrelativistic quantum mechanics on the discretized half-line, 

'— constructing a one-parameter family of Hamiltonians that are analogous to the 
. Robin family of boundary conditions in continuum half-line quantum mechan- 

ics. For classically singular Hamiltonians, the construction provides a singularity 
qq \ avoidance mechanism that has qualitative similarities with singularity avoidance 

encountered in loop quantum gravity. Applications include the free particle, 
the attractive Coulomb potential, the scale invariant potential and a black hole 
described in terms of the Einstein- Rosen wormhole throat. The spectrum is an- 
psj ! alyzed by analytic and numerical techniques. In the continuum limit, the full 

Robin family of boundary conditions can be recovered via a suitable fine-tuning 
but the Dirichlet-type boundary condition emerges as generic. 
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1 Introduction 



It is often useful and sometimes necessary to study quantum mechanics on the 
half-line [0, oo). The boundary, which has been taken to be at the origin without 
loss of generality, can represent the location of a singularity in the potential or 
an infinite potential barrier (impenetrable wall) in the presence of an otherwise 
non-singular potential. The form of the Hamiltonian determines what conditions 
on the wave functions at the boundary lead to a quantization with unitary time 
evolution, if indeed such quantizations exist. Textbook accounts can be found in 
[U [2] and accessible reviews in [5] . 

If the Hamiltonian is essentially self-adjoint, no boundary conditions are re- 
quired. It may however happen that the Hamiltonian has a family of self-adjoint 
extensions, each specified by a boundary condition. As an example, consider a 
particle in a smooth potential. If the potential is sufficiently well-behaved at 
oc y oo so that no boundary conditions are needed there, the self-adjoint exten- 
sions are parametrized by a parameter L G 1U {oo} ~ S* 1 ~ U(l), such that the 
corresponding boundary condition on the wave function is the Robin boundary 
condition, 



The special case L = gives the Dirichlet boundary condition, ip(0) = 0, and 
the special case L = oo gives the Neumann boundary condition, ip'(0) = 0. All 
values of L make the time evolution unitary, and L may be physically interpreted 
in terms of the phase shift of a wave packet on its reflection at the origin. 

From a purely mathematical point of view all the self-adjoint extensions are 
created equal, none are preferred. In most physical situations, however, Dirichlet 
boundary conditions seem to be considered 'natural'. One argument to this end 
comes from regarding a singular potential or an infinite wall as a mathematical 
idealization of a physical situation in which the potential is regular jl] . Assuming 
that the Hamiltonian for the regular potential is already essentially self-adjoint 
and requires no boundary conditions, the 'physical' boundary conditions in the 
singular potential approximation may then be defined as those obtained when 
the regular potential approaches the singular one. It is a somewhat surprising 
mathematical fact that recovering in this limit anything other than the Dirichlet 
boundary conditions tends to require severe fine-tuning in the way in which the 
limit is taken [I]. 

The context of the foregoing discussion was conventional Schrodinger quan- 
tization. In this paper we investigate a similar issue of boundary conditions in 
polymer quantization, in which the continuous variable x G M + U {0} is replaced 
by a discrete variable and the inner product becomes a discrete sum rather than 
an integral [6] . This quantization scheme arises most naturally in loop quan- 
tum gravity [8] . The term "polymer" stems from that fact that the quantum 
states of the so-called "polymer particle represention" are mathematically ana- 



^(0) +Lip'(0) = . 
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loguous to the polymer-like excitations of quantum geometry [5]. Although we 
retain the loop quantum gravity motivated name, we stress that polymer quanti- 
zation is a viable quantization procedure in its own right, unitarily inequivalent to 
Schrodinger quantization. It has been applied to a variety of systems, including 
the harmonic oscillator [S], the Coulomb potential [9], the scale invariant poten- 
tial [10], and the Schwarzschild spacetime as described in terms of the dynamics 
of the Einstein- Rosen wormhole throat [111 [12]. In the last three of these systems 
the configuration space is the half-line, rather than the full real line, and the anal- 
yses in [PJ EU [UJ EE2] treat the restriction to a half-line by first polymer quantizing 
on the full real line and then considering separately the even and odd sectors, 
mimicking respectively the Neumann and Dirichlet boundary conditions in con- 
tinuum. In the case of the even sector, a potential that is singular at the origin 
needs in addition to be regulated, and this was done by the Thiemann trick [13] . 
Most recently, Dirichlet boundary conditions were applied in the analysis of the 
statistical thermodynamics of a polymer particle in a box |14j . 

In this paper we construct a family of polymer Hamiltonians that is analogous 
to the Robin family of boundary conditions in Schrodinger quantization, and we 
analyze the continuum limit of these polymer theories. The construction does 
not rely on extending the polymer dynamics from half-line to full line, and it 
applies to potentials that are singular at the origin without a need to introduce 
a separate regularization for the potential. We shall show that the continuum 
limit of our family of polymer Hamiltonians does reproduce the Robin family 
of boundary conditions for the Schrodinger Hamiltonian. We shall further find 
that a fine-tuning of the polymer Hamiltonian in this limit is required in order 
to recover Schrodinger boundary conditions other than Dirichlet. The continuum 
Dirichlet boundary conditions hence emerge as generic when approached from 
polymer quantization, just as they emerge as generic when approached via regular 
potentials within Schrodinger quantization. 

As concrete examples, we consider the following one-dimensional systems, 
each of which is motivated by specific physical applications: 

1. The free particle in the presence of an infinitely high potential wall. This 
is the prototype for any discussion of quantum mechanics on the half line. 

2. The attractive Coulomb potential. This describes the spherically symmetric 
sector of the hydrogen atom. 

3. The scale invariant potential. This has applications in physics of polymers 
[T5] and in molecular physics [16], and it also arises in the context of quan- 
tizing spherically symmetric black holes [17j . 

4. The Hamiltonian that describes the dynamics of the Einstein- Rosen worm- 
hole throat [TT1[T2]. 
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The paper is organized as follows: In Section [2] we present the general con- 
struction of our half-line polymer Hamiltonians for systems whose classical Hamil- 
tonian is of the conventional nonrelativistic form, H = p 2 + V(x). Sections |3l 
H] and apply this construction respectively to the free particle, the attractive 
Coulomb potential and the scale invariant potential. Section O addresses the 
Einstein-Rosen wormhole throat dynamics: in this system the classical kinetic 
term has a more complicated form, but we show that the construction of Section 
[2] has nevertheless a natural generalization. Section [7] presents a summary and 
brief concluding remarks. 

We use throughout dimensionless units: for the Coulomb potential the con- 
figuration space coordinate is twice the Rydberg radial coordinate as in [9], and 
in the Einstein- Rosen throat theory the units are the Planck units as in (TTJ [12]. 
Complex conjugation is denoted by overline. 



2 Polymer half-line quantization for nonrelativis- 
tic Hamiltonians 

In this section we present our family of half-line polymer Hamiltonians for sys- 
tems whose classical Hamiltonian consists of kinetic and potential terms of non- 
relativistic standard form. Subsection 12.11 recalls the key features of polymer 
quantization of this system on the full real line [5] . The half-line version is given 
in subsection 12.21 



2.1 Polymer real line 

We consider a system whose classical phase space is 1R 2 = {(x,p)} with the 
Poisson bracket {x,p} = 1 and the Hamiltonian 

H = p 2 + V(x). (2.1) 

We assume the potential V to be defined for all x and real-valued. 
The polymer Hilbert space is spanned by the basis states 

*-<*>={£ (2 - 2) 

with the inner product 

(Vk, Vv) = 8x,x>, (2.3) 

where the object on the right hand side is the Kronecker delta. The position 
operator x acts by pointwise multiplication as 

(xt/i)(x) :— xip(x), (2.4) 
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and the family of translation operators | pi G M} is defined as in Schrodinger 
quantization by 

(U^){x) :=il>(x + ii). (2.5) 

In Schrodinger quantization the action of ( 12. 5 p is weakly continuous in /i, 
and the usual momentum operator p := — id x can be defined in terms of the 
translation operator by = e 1 ^ . In polymer quantization the action of is 
not weakly continuous in /i and a similar definition of a momentum operator is 
not available. Instead, we now take fi to be a fixed positive constant, understood 
as a fundamental length scale in the polymer quantum theory, and define the 
momentum operator and its square by 

P:=^(f>,-^), (2.6a) 

p2-.= -^(2-U p -Ul), ( 2 - 6b ) 

where fi := 2fi. The polymer Hamiltonian is defined to be 

H:=f + V, (2.7) 
where V acts by pointwise multiplication, 

(Vi/j)(x) := V{x)^{x), (2.8) 

and 



f = a-Ul). (2.9) 



H is clearly symmetric. Its action decomposes the polymer Hilbert space into 
a continuum of separable superselection sectors, where each sector consists of wave 
functions supported on the regular /i-spaced lattice {A + mfi m 6 Z} and the 
parameter A e [0,fi) parametrizes the sectors. A study of the deficiency indices 
[1] shows that H has at least one self-adjoint extension in each superselection 
sector. In particular, the Kato-Rellich theorem [1] can be applied as in [10] to 
show that H is essentially self-adjoint in every superselection sector in which the 
function V is bounded. 

It is useful to write the Hilbert space of each superselection sector in terms 
of sequences. For concreteness, consider the superselection sector A = 0. An 
orthonormal basis is given by {ip m ft | m G Z}. Writing ip = J2 m c m ip m fi, the 

inner product reads (tjj^, ip^~) = c m W c^', and the action of H reads 

H [y^Cm^mp J = ^ (— "^f + V(mfi)CrnJ ^mp, ■ (2.10) 
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The Hilbert space hence consists of the two-sided square summable sequences 

oo 

111= — oo ' 



c := (c m )^ = , and the action of H reads 



(Hc) m = 2Cm ^ Cm ~ X + V(mfi)c m . (2.11) 
2.2 Polymer half-line 

We wish to modify the polymer Hamiltonian (12. lip into a Hamiltonian in the 
Hilbert space of one-sided square-summable sequences. 

We take the one-sided sequences to be of the form c := (c m )^ =1 and the inner 
product to read (d,c) = Ylm=i^ mCm - We define the one-parameter family of 
modified Hamiltonians \H a \ a E Mj by 



2c m c m _|_i c 



m—l 



r,2 



+ V(m/j,)c m for m > 1; 



{H a c) m :={ * (2.12) 
1 (2-a)ci-pa . . 
— h k (w)ci for m = 1. 

Each if a is symmetric, and it can be verified as with the two-sided sequences 
that each H a has at least one self-adjoint extension. If V = 0, it can be explicitly 
verified that there are no nonvanishing normalizable solutions to H a c = ±ic, and 
in this case H Q is hence essentially self-adjoint. It follows by the Kato-Rellich 
theorem [1] that H a is essentially self-adjoint whenever the set (V(mp) \ m = 
1,2,...} is bounded. 

The motivation for the definition of (H a c) by fl 2 . 1 2 [) is that (H a c) agrees 
with [He) (12. lip for m > 1 if we envisage there to be a fictitious lattice point 
at m = such that c = ac\. This means that a e K can be regarded as a lattice 
counterpart of the continuum theory self-adjoint extension parameter L. The 
case a = is analogous to continuum Dirichlet, with the continuum boundary at 
m = 0. The cases a = ±1 are respectively analogous to the continuum Neumann 
and continuum Dirichlet, with the continuum boundary half-way between m = 
and m — l. 

Two comments are in order. First, given the definition (12.121) . the coefficient 
of d\ c\ in the inner product (d, c) = Ylm=i c m cannot be changed without 
losing symmetricity of H a . Second, although the definition of H a is motivated 
by a fictitious lattice point at m — 0, the actual definition involves V only at 
lattice points x = mfi with m > 1. H a is hence well defined even if the classical 
potential is singular at x — 0. 



3 Free particle 

As the first example we consider the free particle, V = 0. We begin by recalling 
relevant facts from Schrodinger quantization on the half-line [H El E] and from 



5 



polymer quantization on the full real line [5]. We then analyze H a (I2.12p and its 
continuum limit. 

3.1 Schrodinger half-line 

In Schrodinger quantization, the Hamiltonian operator of a free particle on the 
half-line is Hf rec := —d 2 /dx 2 , acting in the Hilbert space L 2 (1R+, dx). Hf rcc is sym- 
metric and has a U(l) ~ S l of self-adjoint extensions, specified by the boundary 
condition (11. ip . 

The spectrum of each extension contains the positive continuum. The exten- 
sions with < L < oo have in addition a discrete ground state of energy —L~ 2 . 
The spectral decomposition of the identity in terms of the eigenfunctions can be 
found in [18]. 

3.2 Polymer real line 

Polymer quantization on the full real line can be analyzed by interpreting the 
two-sided square summable sequence c := (c m )^ = _ oc as the coefficients in the 

Fourier expansion x(<p) = {^Y 1 ^ 2 Y^ m c m^ mip of the function \ L 2 (S 1 ), so 
that in this realization (Hx)(<p) = 2/i~ 2 (1 — cos<^)x(y9). It is immediate that H 
is essentially self-adjoint and its spectrum is the continuum (0, 4/i -2 ). 

3.3 Polymer half-line 

The half-line polymer Hamiltonian is given by H a ( 12.12ft with V — 0. As noted 
above, H a is essentially self-adjoint. 

3.3.1 Spectrum 

We wish to find the spectrum. 

Let E G M denote the energy eigenvalue, which may be a proper eigenvalue 
in the discrete spectrum or an improper eigenvalue in the continuous spectrum. 
The eigenvalue equation reads 



If E 7^ and E ^ 4/x 2 , the solution to (13.11) as a difference equation is 



H a c = Ec . 



(3.1) 



in 



aA™ + bA 



(3.2) 



where a and b are constants, 



A± := (1 - \fE) ± ^/(1-l^E) 2 -1 , 



(3.3) 
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and evaluating (13. ip at the lattice point m = 1 by (12.121) shows that (13.21) must 
in addition satisfy 




c 2 = (2 - a - Ji 2 E) a , (3.4) 

so that for each a there is only one linearly independent solution. If E = 
or E = 4:fi~ 2 , we proceed similarly, finding that the only linearly independent 
solution to (13.11) as a difference equation is 

a)m for E = 0; 

(3.5) 

- (1 + a)m] for E = 4/i" 2 . 

The spectrum is now found by considering the normalizability properties of these 
solutions. 

Suppose first that < E < Afl~ 2 . We parametrize E as 

tfE = 4sin 2 (fl/2) , (3.6) 

where < 9 < tt. From (13.31) we then have A± = e ±te , and from (13.41) it follows 
that the solution reads 

c m = sm(m9 + 5) , (3.7) 
where 5 is determined in terms of a and 9 from 

cot 5 = — —-cote, (3.8) 
a sin 9 

understood for a = in the limiting sense 5 = 0. Note that 5 is unique up to an 
additive integer multiple of tt. Note also that in the special cases a = 0, a = 1 
and a = —1, we have respectively 5 = 0, 8 = {ji — 9)/2 and 8 = —9/2. As the 
solutions are oscillatory at m — > oo, they are lattice analogues of plane waves, 
belonging to the continuous spectrum. 

Suppose then that E < 0. We now have A + > 1 and < A_ < 1, and 
avoiding an exponential divergence at m — > oo in (13. 2p requires a = 0. Matching 
to ( 13. 4 p is possible iff a > 1, and E is then uniquely determined in terms of a by 

E= (2-a-a" 1 )/i~ 2 . (3.9) 

The solution is the normalizable, proper eigenstate c m = a~ m . 

Suppose next that E > 4/i -2 . A similar analysis shows that a solution exists 
iff a < —1, E is uniquely determined in terms of a by ( I3.9p . and the solution is 
the normalizable, proper eigenstate c m = a~ m . 

Suppose finally that E = or E = 4/i -2 . As the solutions (I3.5P are not 
normalizable, these special values of E are not proper eigenvalues, and as points 
of measure zero they do not contribute to the spectrum. 
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We summarize. For any a, the spectrum contains the band (0,4/i -2 ) of im- 
proper, Dirac-delta-normalizable eigenstates, given by (13. 7p with (13. 6p and (13. 8p . 
For each a such that \a\ > 1, there exists in addition a unique proper eigenstate, 
given by c m = a~ m and having the eigenvalue (13. 9p . For a > 1 the proper eigen- 
state is a ground state at E < 0, below the continuous band, and for a < — 1 it 
is a highly excited state at E > 4:Ji~ 2 , above the continuous band. 

We note that the special case a = 1 of our free particle Hamiltonian is obtained 
from the family of discrete Hamiltonians (2.12) in [19] with the choices B(x) = 
D(x) = 1 for x > and B(0) = 1, where the discrete configuration variable x 
of [IS] is such that x = corresponds to our m = 1. As discussed above, the 
spectrum with this value of a consists of the continuum band (0, 4/2 -2 ) and does 
not contain a bound state. This exemplifies the observation made at the end of 
subsection 12.21 that even though we invoked a fictitious m = lattice point to 
motivate the definition of the Hamiltonian, the actual definition does not refer to 
this lattice point and the value of the wave function at the fictitious lattice point 
never enters the theory: in the formulation of [I2| the a = 1 theory is indeed 
constructed without a motivation via fictitious lattice points. 



3.3.2 Continuum limit 



We now show that the full one-parameter family of continuum Schrodinger quan- 
tizations can be recovered in the continuum limit. 

The key for taking the limit is to make a dependent on Jx in a suitable way. 
We introduce a parameter Lp e R U {oo}, which is considered independent of fi, 
and we choose 

'0 forL P = 0; 



a 



for Li 



oo; 



(3.10) 



otherwise, 



. 1 " ifl/Lp) 

which is well defined for sufficiently small Jx when Lp is in the interval (0, oo), and 
for all Jl when Lp takes other values. When interpreted in terms of the fictitious 
lattice point at m = for which Cq = aci, (I3.10p amounts to 



Co + Lp - 1 — — = for Lp ^ oo; 



/' 

(ci - c ) 

ft 



(3.11) 



for Lp = oo. 



Comparison of (ll.ip and (13.111) suggests that the polymer theory should reduce 
to the continuum theory with L = Lp as Jl — > 0. We shall verify that it does. 

Let first < Lp < oo and focus on the ground state. By (13. 9p . the ground 
state energy is given by E = —Lp 2 + O(Jl), and the ground state wave function 
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can be written in terms of the discrete distance coordinate x = m\x as 

tip) ■= = (l - " = e- xjLp + 0(fi), (3.12) 

where in the last expression Lp and x are considered independent of p. In the 
limit p — > 0, we therefore recover the ground state energy and the ground state 
wave function of the continuum theory with L = Lp. 

Let then Lp be arbitrary and consider energies in the continuous band. For 
fixed E, (ESD and (E2D give 

{0 for L P = 0; 

§7T + 0(/2) forL P = oo; (3.13) 
— arctan^LpA/E j + 0{p) otherwise, 

and 

ip(x) := c x/A = sin((9/p)x + 5) = sin(VEx + 5) + 0{fi 3 ). (3.14) 

In the limit Jx — > 0, we therefore recover the continuum theory wave functions 
with L = Lp. 

We emphasize that the fine-tuning of a as a function of p, in ( 13 . 1 j) is essential. 
If the limit Jx — > is instead taken with fixed a, the choice a = 1 (corresponding 
to Lp = oo in (13.101) ) gives the Neumann continuum theory, L = oo, while 
any other choice for a (corresponding to Lp — >■ in (13.10jl ) gives the Dirichlet 
continuum theory, L = 0. 



4 Coulomb potential 

As the next example, we consider the attractive Coulomb potential, V(x) = —1/x. 
A new feature is that this potential is singular at x — 0. 



4.1 Schrodinger half-line 

In Schrodinger quantization [20], the Hamiltonian has a one-parameter family of 
self-adjoint extensions, specified by a boundary condition at x = 0. The boundary 
condition does not take the Robin form (11. ip but reads instead 

(L - (1) - (2)) ^(0) + lim f + M \ i nx ) = 0) (4.1) 

x^>0 \ X J 

where ip is the wave function, \1/ is the digamma function [21] and L G K. U {oo} 
is the parameter. The second term in (14.11) can be understood as a logarithm- 
corrected first derivative, where the correction must be included because of the 
singularity of the potential at x — > 0. 
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The spectrum consists of the positive continuum and a countable number of 
proper negative eigenvalues. Writing s := 1/v— 4E for E < 0, the eigenvalues 
are given by the solutions to 

L = *(l-s)-lns + — . (4.2) 

The solutions to (14. 2 p with L = oo are s = 1, 2, 3, . . ., coming from the simple 
poles of ^ at non-positive integers. This extension is analogous to the Dirichlet 
boundary condition for regular potentials, and it arises via a limiting prescription 
from regularized potentials that make the three-dimensional attractive Coulomb 
Hamiltonian essentially self-adjoint [2D]. This is the extension that yields the 
textbook spectrum, realized in nature. 

We note in passing that if physical dimensions are restored, L has an inter- 
pretation as an inverse length scale. Using L, rather than the associated length 
scale 1/L [2D], is convenient for us because the spectrum changes smoothly across 
L = 0. The only discontinuous change in the spectrum in terms of L occurs when 
L — > oo through positive values: in this limit the ground state descends to — oo 
and disappears. 



4.2 Polymer half-line 

The half-line polymer Hamiltonian is given by H a (I2.12p with V{mji) = — 1/ (m/x). 
As noted in subsection 12.2} the boundedness of V over the lattice points shows 
that H a is essentially self-adjoint. 

We focus on a numerical study of the proper eigenvalues in the limit ]x — > 0. 
The numerical scheme is as described in [9J except that the boundary condition 
to shoot for is now c = olc\. 

The key is again to make a depend on Ji in a suitable way. Motivated by the 
continuum boundary condition (14. ip . we introduce a parameter Lp G R U {oo}, 
and we choose 

{0 for Lp = oo; 
1 . . (4-3) 
-= otherwise, 
l-/2(L P -*(l)-*(2)+ln/z) 

which is well defined for sufficiently small Ji with any fixed value of Lp. Numerical 
results for the lowest three eigenvalues in terms of Lp are shown in Table [TJ 
Numerical accuracy allows us to probe values of ft down to 10~ 5 . 

For Lp = oo, we find that the polymer eigenenergies converge to those of the 
Schrodinger theory with L = oo. For other values of Lp, we find that the polymer 
eigenenergies converge to those of a Schrodinger theory with L — Lp w 0.423. 
While we do not have an analytic explanation for this shift in L — Lp, we note 
that within our numerical accuracy the shift coincides with the constant ^/(2) w 
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2.100712343 


2.4607015522 


2.54141508 


2.63445759 


2.90674777 


34 


1.042 x 10" 


10 


P = 10- 5 


2.100712544 


2.4606208222 


2.54132730 


2.63437572 


2.90673263 


34 


1.042 x 10" 


12 


Schrodinger 


2.09672 


2.4193 


2.4987 


2.5970 


2.9030 


3 






L p — L c ft 
p = ICT 3 


0.420509240 


0.4296974160 


0.4306588703 


0.431510579 


0.44284175 








M = 10" 4 


0.422787147 


0.4237079042 


0.4237880518 


0.423888137 


0.42479883 








p = It)" 5 


0.422807636 


0.4228997124 


0.4229077629 


0.422917713 


0.42300750 









Table 1: The table shows numerical results for the lowest three eigenenergies in 
the Coulomb potential, in the polymer theory as a function of the parameter Lp 
and the scale fi, and in the Schrodinger theory as a function of the parameter L. 
The shown quantities Sq, Si and S2 parametrize the eigenenergies E , E\ and E 2 
by s = 1 / yJ—AE. The last three rows for each Sj show the values of Lp — L e g such 
that the Schrodinger eigenvalue with L = L e g equals the corresponding polymer 
eigenenergy. Note the convergence of Lp — L e g to 0.423 with decreasing ft, within 
the numerical accuracy. 
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0.4227843351 that appears in (14. ip and (14. 3p . Examining whether this coincidence 
continues to hold beyond three decimal places would require numerical work 
beyond ft = 10~ 5 . 

If the limit p, — y is taken with any fixed value of a, the relation f)4.3p 
shows that the numerical results in Table [T] are consistent with convergence to 
the continuum theory with L = oo. 



5 Scale invariant potential 

We next consider the scale invariant potential, V(x) = — X/x 2 , where A is a 
real- valued constant. Although we work in dimensionless variables, we note that 
A remains dimensionless even when physical dimensions are restored, and A is 
hence a pure number whose value is significant regardless any unit choices. That 
the coupling constant is dimensionless is the special property of a scale invariant 
potential. 



5.1 Schrodinger half-line 

Schrodinger quantization of the system is reviewed in [TO]. For A < —3/4 the 
Hamiltonian is essentially self-adjoint and the spectrum consists of the positive 
continuum. For A > —3/4 the self-adjoint extensions are specified by a parameter 
that takes values in U(l) ~ S 1 , and the spectrum of each extension contains the 
positive continuum but discrete negative eigenenergies can also exist. We shall 
recall here relevant facts about the cases in which discrete eigenenergies do exist. 

For A > 1/4, each extension has a countable tower of proper eigenenergies, 
given by 



E n = - exp 



2( 7 -7rn)/VA- (1/4)1, n G Z, (5.1) 



where 7 G [0, 7r) is the parameter specifying the extension. The discrete spectrum 
is unbounded below, E n — > —00 asn-y —00, and it accumulates to from below, 
E n — )■ 0_ as n — > 00. 

For A = 1/4, one of the extensions has no discrete eigenenergies but the rest 
have exactly one negative eigenenergy each. 

For —3/4 < A < 1/4, there is an open interval of extensions that have ex- 
actly one negative eigenenergy each. The remaining extensions have no discrete 
eigenenergies. 



5.2 Polymer half-line 

The half-line polymer Hamiltonian is given by H a (I2.12p with V(m/2) = — A/(m/i) 2 . 
The boundedness of V over the lattice points again implies that H a is essentially 
self-adjoint. 
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We again focus on a numerical study of the discrete eigenvalues. Because of 
the scale invariance of the potential, the polymer eigenvalues depend on p, only 
through the overall factor fir 2 . The question to be examined hence is how the 
combination fi 2 E depends on A. 

Because of the strong singularity of the continuum potential at x — > 0, it is not 
clear how the Schrodinger theory self-adjointness boundary conditions [10] might 
motivate a value of a in our polymer Hamiltonian. We hence take a = tanx, 
where x [0,vr), allowing x — 7r /2 to be considered as a limiting case. The 
numerical implementation is as in [TU] except that the boundary condition to 
shoot for is 



For A > 2, we find a tower of negative eigenenergies E < E 1 < E 2 < ■ ■ ■ < 0. 
We have followed the tower up to Eq but slowness of the numerics has not enabled 
us to examine whether the tower terminates. Starting at x — and increasing x, 
E first migrates downwards and then disappears at x ~ 7r /2; slowness of the 
numerics has not enabled us to examine whether this disappearance happens by 
descent to — oo or by some other mechanism. The excited states E n , n > 0, 
migrate down smoothly, meeting at x ~ > 71 with E n _i at % = 0. Graphs for E 3 
in terms of x are shown in Figured] for A = 4 and A = 8. 

When A decreases towards 0, the numerics becomes slow. Outside the interval 
7r/4 < x < tt/2 we have found no evidence that any eigenstates would survive. 
Within the interval ir/4 < x < 7r /2, however, comparison with the free particle 
polymer theory of Section [3] suggests that the ground state should persist as 
A —> and tend to that of the corresponding A = theory, given by 



and we find numerical evidence that this indeed happens. A sample plot is shown 
in Figure |2j 

If the lowest few states are excluded, the eigenergies for A > 2 are a good fit to 
the continuum eigenenergy formula (15. ip , as shown in Figure [3j The continuum 
parameter 7 is an increasing function of Xi an d for large A it is slowly varying 
except when x is close to tt/2. 

6 Einstein-Rosen wormhole throat 

As the final example, we study the interior dynamics of the eternal Schwarzschild 
black hole, adopting as the configuration variable the area of the Einstein-Rosen 
wormhole throat and as the time parameter the proper time of a comoving ob- 
server at the throat [221 [23| 124] . 



c cos x — C\ sin x = 0. 



(5.2) 




(5.3) 
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S J 




Figure 1: Scale invariant potential. S3 := — ln(— fi 2 Es) is shown as a function of 
X for A = 4 (upper curve, red) and for A = 8 (lower curve, blue). 
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Figure 2: Scale invariant potential. Values of s n := — ln(— fi 2 E n ) are shown as 
a function of A for \ = 7r/3, for n = (low, black), n = 1 (middle, blue) and 
n = 2 (top, red). As A — > 0, s approaches the value of the A = polymer theory, 
— ln(— 2 + 4/V3) ~ 1.11731, and the excited states disappear. 
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2 3 4 5 6 

n 

10- 
8- 

s_n 6 _ 

4- 
2- 

T ■ 1 ■ 1 ■ 1 ■ 1 

2 3 4 5 6 

n 

Figure 3: Scale invariant potential. Values of s n := —\n(—fi 2 E n ), 2 < n < 6, 
are shown for different values of Xi with a linear regression fit to the continuum 
formula ( 15. 1ft . Upper diagram shows A = 4 and the lines from top to bottom are 
with x — (red), x — 0.557T (blue), x — 0.607T (brown) and x = 0.657T (green). 
Lower diagram shows A = 8 and the lines from top to bottom are with x — 
(red), x = 0.557T (blue) and x = 0.567T (brown). 
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6.1 Schrodinger half-line 

Following pX], we take the classical phase space to be (<f), II), where the positive- 
valued configuration coordinate is the square of the wormhole throat area-radius 
and II is the real- valued conjugate momentum. The classical Hamiltonian is given 
by 

H=l^ 2 (4U 2 + l) . (6.1) 

We take the wave functions in Schrodinger quantization to be square integrable 
functions of in the measure d<p. For reasons that will emerge in the polymer 
theory, we consider the quantum Hamiltonian with the symmetric ordering 

^=H- 4 ^ +#/2 )- (6 ' 2) 

It is shown in [TTJ [21] that the self-adjoint extensions of -f/schr are labelled by the 
parameter 9 G [0, ir) such that 

^,(0) = [1 + O(0 3/2 )] sin 9 + 1/2 [1 + O(0 3/2 )] cos 6 , (6.3) 

or equivalently 

^(0) cos 9 - ( lim m -? {0) ) sin 9 = , (6.4) 

so that 9 = is analogous to Dirichlet and 9 = tt/2 is analogous to Neumann. The 
spectrum of each self-adjoint extension is discrete and bounded below. Further 
analytic and numerical results on the spectrum are given in [HI [21] • 

6.2 Polymer real line 

As the kinetic term of the classical Hamiltonian (16. ip does not have the conven- 
tional form of (12.11) . the polymer quantization prescription of Section [2] must be 
suitably generalized. We now recall the results of [11] for the situation where <fi 
is polymer quantized on the full real line. 

We polymerize the pair (0,11) as the pair (x,p) in subsection 12.11 letting 
take all real values, defining the basis states 

with the inner product 

(^>,Vy) = <W'> ( 6 - 6 ) 
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and the operators 

(6.7a) 

(^)(0):=</# + /i), (6.7b) 
ft-^ M -ej), (6.7c) 

where the positive constant /i is again the fundamental polymer length scale. The 
polymer counterpart of the Schrodinger Hamiltonian f)6.8p is taken to be 

H=-Uil^ 2 Ti + ^ 2 ). (6.8) 



2 

where the operator 1 / 2 is defined by 

(0 1 / 2 ^)(0):=|0| 1/2 ^(0). (6.9) 

Specialising to the superselection sector in which is confined to the lattice 
[mji meZ} with p, := 2/x, and writing ip = J^m^^mp,, the action of H on 
the two-sided square-summable sequence c := (c m )°? reads 



(*<=)» = ^[(h + ir + h-il^ + W* 

I i I 1/2 I i I 1/2 

\ m ~T 2 I c m+l ~ | m _ 2 I C m-1 



(6.10) 



It follows that the polymer Hilbert space breaks further into the even superse- 
lection sector, in which c m = c_ m , and the odd superselection sector, in which 
c m = — c_ m . In the limit Jx — > 0, it is found numerically [11] that the even (re- 
spectively odd) sector converges to the Schrodinger theory with 9 = n/2 (6 = 0). 
This convergence is consistent with what one would expect by inspection of the 
Schrodinger theory boundary condition ( 16. 3p . 



6.3 Polymer half- line 

We wish to modify the polymer Hamiltonian (I6.10p to act on one-sided square- 
summable sequences, of the form c := (c m )^ =1 and with the inner product 
(d,c) = J2m=idm,c m - We again choose the one-parameter family of modified 
Hamiltonians \H a \ a 6 Rj to be defined by adding at m = a fictitious lattice 
point with with cq = aci, so that 



(He) 



3/2 



(m + l) 1/2 +(m-|) 1/2 + \fm 1 / 2 ) 



- (m + |) 1/2 c m+1 - (m - l) 1/z c m _i for m > 1, (6.11a) 



ni/2 



3/2 



nV2 



'3\l/2 



C 2 



(6.11b) 
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6 P or 





tt/4 


tt/2 


3?r/5 


4tt/5 


ft=l 






1.23093 


1.11371 


0.904625 


0.592200 




ft = 10- 


1 




1.18178 


1.04988 


0.785426 


0.423878 




ft = 10- 


2 




1.16602 


1.03780 


0.770035 


0.388607 




/i = io- 


3 




1.16113 


1.03478 


0.768373 


0.380709 




/x = 10- 


4 




1.15960 


1.03390 


0.768203 


0.378624 




Schrodin^ 


cr 


1.15890 


1 


0.768184 


0.535489 


-3.14921 







0.19527 


1/2 


0.636933 






9 P or 





tt/4 


tt/2 


3tt/5 


4tt/5 


M = l 






1.87343 


1.80220 


1.72038 


1.65342 


1.39014 


ix = 10- 


1 




1.84343 


1.75446 


1.63581 


1.54923 


1.32377 


fi = io- 


2 




1.83184 


1.74589 


1.62456 


1.53274 


1.30281 


/i = io- 


3 




1.82825 


1.74396 


1.62336 


1.53023 


1.29678 


Al = 10- 


4 




1.82713 


1.74343 


1.62324 


1.52975 


1.29496 


Schroding 


er 


1.82661 


1.72401 


1.62322 


1.56006 


1.34989 







0.19526 


1/2 


0.63692 


0.84817 


S 2 


6» P or 





tt/4 


tt/2 


3?r/5 


4tt/5 


ft=l 






2.34359 


2.28992 


2.23644 


2.19574 


2.00416 


ft = io- 


1 




2.32411 


2.25111 


2.16639 


2.11021 


1.95368 


Al = 10- 


2 




2.31424 


2.24373 


2.15637 


2.09634 


1.93652 


/x = 10- 


3 




2.31117 


2.24216 


2.15530 


2.09435 


1.93177 


Al = 10- 


4 




2.31022 


2.24174 


2.15519 


2.09401 


1.93036 


Schrodin^ 


er 


2.30978 


2.22684 


2.15518 


2.11361 


1.97162 







0.19524 


1/2 


0.63691 


0.84817 


E:>, 


P or 





tt/4 


tt/2 


3?r/5 


4tt/5 


ft = l 






2.73280 


2.68919 


2.64913 


2.61958 


2.46137 


ft = io- 


1 




2.72106 


2.65720 


2.58907 


2.54603 


2.42038 


Al = 10- 


2 




2.71224 


2.65049 


2.57977 


2.53351 


2.40498 


Al = 10- 


3 




2.70948 


2.64911 


2.57877 


2.53177 


2.40087 


Al = 10- 


4 




2.70862 


2.64875 


2.57867 


2.53148 


2.39961 


Schroding 


er 


2.70822 


2.63625 


2.57866 


2.54653 


2.43449 


e g/ 7T 





0.19526 


1/2 


0.63691 


0.84816 



Table 2: The table shows the four lowest eigenenergies for the Einstein- Rosen 
wormhole throat, in the polymer theory as a function of the parameter 9p and 
the scale p,, and in the Schrodinger theory as a function of the parameter 9. For 
each dp, the polymer eigenenergies converge to those of the Schrodinger theory 
with 9 = 9 c g as shown. 
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A(e) 




Figure 4: A(9) := (dp — O c g) is shown as a function of Op for the Einstein-Rosen 
wormhole throat. 



Each H a is symmetric, and the coefficient of d\ C\ in the inner product cannot 
be changed without losing this symmetricity. The growth of the potential term 
suggests, by comparison with the solutions to the continuum eigenvalue equation, 
that each H a is essentially self-adjoint; we have however not attempted to examine 
this rigorously. 

We wish to examine the bound states in the limit Ji — > 0. As in the previous 
cases, the key is to make a depend on ft in a suitable way. Motivated by the 
continuum boundary condition (16. 4p . we introduce a parameter 9 P G [0, 7r), and 
we choose 



which is well defined for sufficiently small ft with any fixed value of 9p. Numerical 
results for the four lowest eigenvalues in terms of Op are shown in Table [2j 

For Op = and Op = tt/2, the polymer eigenenergies converge to those of 
the Schrodinger theory with respectively = and = n/2: this reproduces 
the results found in the even and odd sectors of the real line polymer theory 
in pX]. For other values of Op the polymer eigenenergies converge to those of 
the Schrodinger theory with = c g, which differs from Op as shown in Figure HI 
Note that when Op runs over the full interval [0, 7r), so does c q. 

If the limit p, — y is taken with a = 1, it is seen from (16.121) and Table [2] 
that the polymer theory converges to the continuum theory with = tt/2, and if 
the limit is taken with any other fixed value of a, the polymer theory converges 
to the continuum theory with = 0. The continuum Dirichlet-type boundary 
condition hence again emerges as generic. 



sin Op 



(6.12) 



a = 



y/JIcosOp + sin^p ' 
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7 Conclusions 



We have constructed a one-parameter family of polymer quantization Hamil- 
tonians on the half-line. This family is mathematically analogous to the one- 
parameter family of Robin boundary conditions in Schrodinger quantization on 
the half-line. For the free particle, the attractive Coulomb potential and the 
Einstein- Rosen wormhole throat, we found that the full family of continuum 
Robin boundary conditions can be recovered in the continuum limit provided the 
polymer parameter is suitably fine-tuned, while without such a fine-tuning the 
continuum limit yields a Dirichlet-type continuum boundary condition. For the 
scale invariant potential the spectrum depends on the polymer scale only by an 
overall scaling, and a proper notion of a continuum limit does hence not exist, but 
even in this case the spectral properties revealed a close correspondence between 
the polymer theory parameter and the continuum theory boundary condition. 

A notable feature of our half-line Hamiltonian is that it is well defined even 
when the classical potential function is singular at the origin, without any need to 
regularize the singularity. Mathematically, this feature arises because even though 
a fictitious lattice point at the origin was invoked to motivate the definition of 
the Hamiltonian, the actual value of the wave function at the origin never enters 
the theory, neither through the kinetic term nor through the potential term in 
the Hamiltonian. Physically, this feature is relevant to the issue of singularity 
resolution in polymer quantum gravity [25l [26] and may also have some impact 
in situations where polymer quantization is adopted gulator of ultraviolet 

divergences in quantum field theory [27]. 

One qualitative difference between the one-parameter family of continuum 
half-line Hamiltonians and our one-parameter family \H a } of polymer half-line 
Hamiltonians is that the continuum parameter takes values in U(l) ~ S 1 ~ 
R U {oo} but our polymer parameter a takes values in R. While the definition 
f)2.12p of H a does as such not have a well-defined limit as a — > ±oo, one could ask 
whether the (generalized) eigenstates of H a might have an a — > ±oo limit that is 
sufficiently regular for the family {H a } to be completed into an S 1 . Heuristically, 

thinking of H a in terms of the fictitious lattice point at m = such that cq = ac\, 
the a — > ±oo limit should mean taking c\ = 0, which in turn should be equivalent 
to Hq on a lattice that is shifted to the right by one step. For the free particle, we 
can verify from the analytic solution of Section [3] that this heuristics is correct: 
as a — > ±oo, ( 13. 9 p shows that the bound state disappears, while (13 .8p shows that 
states in the continuous spectrum have the limit 5 — > —6, which gives exactly 
the expected shift in m in (13. 7p . It would be interesting to examine whether the 
heuristics holds also for nonconstant potentials. 

In all our examples, the parameter in the half-line polymer Hamiltonian af- 
fects the qualitative properties of the spectrum in a way that mimics closely the 
effects of the self-adjointness parameter in the continuum half-line theory. For 
example, for the free particle we saw that a single normalizable state can be made 
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to appear, and the energy of this state can be tuned to any value outside the con- 
tinuous spectrum. This is just as in the continuum theory, the only difference 
being that the polymer continuous spectrum is an interval while the continuum 
continuous spectrum is half-infinite. As another example, we saw instances where 
adjusting the polymer parameter makes the lowest energy eigenvalue disappear 
while the higher eigenvalues each migrate downwards by one notch: similar be- 
haviour occurs in the continuum theory. We were able to give a fully analytic 
description of these phenomena only for the free particle, but our numerical re- 
sults suggest that it would be interesting to pursue an analytic description also 
for nonconstant potentials. 

The results of our case studies are mutually compatible, and compatible with 
the continuum Schrodinger theory, to an extent that may exemplify a generic 
behaviour of the half-line polymer theory as a function of the parameter in the 
Hamiltonian. To examine this further, it would of course be of interest to extend 
the case case studies to wider classes of classical potentials, controlling both the 
analytic and numeric aspects. For example, one might want to consider the 
harmonic oscillator, for which the Schrodinger problem is readily analytically 
solvable in terms of parabolic cylinder functions [21]. We expect the half-line 
harmonic oscillator polymer theory to be related to the Schrodinger theory in a 
way qualitatively similar to what we found for the wormhole throat; however, we 
have not been able to develop a sufficient numerical control of the polymer theory 
eigenvalues to examine this question quantitatively. 
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